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Abstract. Let ft be a symmetric cone. In this note, we introduce 
Hilbert's projective metric on Q in terms of Jordan algebras and 
we apply it to prove that given a linear transformation g such that 
g(Q) C fl and a real number p, \p\ > 1, then there exists a unique 
element x S f2 satisfying g{x) = x p . 



1. Introduction 

If A is a positive linear transformation on M. n , then the Perron The- 
orem says that there exists xq G M™ \ {0} such that for all x G R™, 
A n x converges in direction to Xq, i.e. ppfj| — ¥ yf^if- ^ ^ as been shown 
by Birkhoff jlj that the Perron Theorem can be considered as a spe- 
cial case of the Banach contraction theorem. Birkhoff 's approach is to 
consider a projective metric in a cone of positive elements in a Banach 
space. This metric was first introduced by Hilbert [S] when studying 
hyperbolic geometry. We begin with the definition of this metric in a 
general setting. 

Let V be is a real Banach space and C be a closed convex pointed 
cone, where pointed means that C PI — C = {0}. Denote the corre- 
sponding interior by Q. The relation =^ is defined on V by saying that 
x =^ y if and only if y — x G C. Hence (V, =4) is a partially ordered 
linear space and it is Archimedean, that is, if ny y for all n G N*, 
then y ^ 0. 

For x G V and y G f2 we let 

M(x,y) := inf{A | x ^ Ay}, 

and 

m(x,y) := sup{/i | /!?/ ^ x}. 
Hilbert's projective metric is defined on Q by 

M(x,y) 



(1.1) d(x,y) = log 



m(x, y) 
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In the case of M™, Hilbert's projective metric is d(x,y) = log . If 

Vi 

where x = (x\, . . . , x n ) and y = (yi, . . . , y n ) are two vectors of IR™ . 

The Hilbert projective metric may be applied to variety of problems 
involving positive matrices and positive integral operators. For example 
one can use it to solve some Volterra equations. It is also particularly 
useful in proving the existence of the fixed point for positive operators 
defined in a Banach space. In this way, it has been shown by Bushell 
j3] that Hilbert projective metric may be applied to prove that, if T a 
real nonsingular r x r matrix, then there exists a unique real positive 
definite symmetric r x r matrix A such that 

(1.2) T'AT = A 2 . 

Notice that if T is neither symmetric nor orthogonal the existence and 
the uniqueness of A is not an elementary problem, even if r = 3. 

In this note we will formulate the Hilbert projective metric on sym- 
metric cones in a way most convenient for our purpose using Jordan 
algebra theory and extend Bushell's Theorem to this class of convex 
cones. 

2. Preparatory definitions 

We summarize here a number of basic definitions and results con- 
cerning symmetric cones. For more detailed discussion, see Faraut and 
Koranyi [5]. 

Let Q be an open convex cone in a Euclidean vector space V of 
dimension n. Let G(Q) be the group of linear automorphisms of Q 

G(Q) = {ge GL(V) \ g{Q) = Q}. 

Then Q is said to be homogeneous if G(Q) acts on it transitively. If 
Q is pointed, then Q is said to be symmetric if it is homogeneous and 
self-dual. 

A Euclidean Jordan algebra is a Euclidean vector space V equipped 
with a bilinear product such that 

xy = yx 

x(x 2 y) = x 2 (xy) 

(xy\z) = (y\xz). 

It is shown that the interior Q of the set of squares in V is a symmetric 
cone, and every symmetric cone is given in this way. 
We define the (left) multiplication L by L(x)y = xy and the so-called 
quadratic representation P by P(x) = 2L 2 (x) — L(x 2 ). For any x G V, 
the endomorphisms L(x) and P(x) are self-adjoint. 
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Example 2.1. The space Sym(r, M) of real symmetric r xr matrices is 
a Jordan algebra for the product x o y = \{xy + yx) and it is Euclidean 
for the scalar product (x\y) = Tr(xy). The associated symmetric cone 
is the set Qsym of positive definite symmetric r xr matrices and G(f2s ym ) 
is the linear group GL(r, M). In this case, P(x)y = xyx. 

The Euclidean Jordan algebra has and identity element denoted by 
e. Let K(Q) := {g G G(Q) \ g(e) = e} = , then K(Q) is a maximal 
compact subgroup of G(Q). Let r be the rank of V. A Jordan frame 
{ci, . . . , ev} of V is a complete system of non-zero orthogonal primitive 
idempotents : 

c\ = Ci, Ci indecomposable, 
CiCj = if % / j, 
ci + . . . + c r = e. 

We denote by J{V) the set of all primitive idempotents of V . Suppose 
V is simple. In other words, there is no non-trivial ideal in V . Then 
each element x in V can be written as 

r 

(2.1) x = kf^XjCj), k e K{Q) , Xj e R, 

3=1 

where K(Q) a is the identity component of K(Q). The determinant is 
defined by det(x) = YYj=i anc l ^ ne t race form by tr(x) = Y7j=i ^j- 
The real numbers Ai, . . . , A r are the eigenvalues of x and the spectral 
norm of x is then defined by 

(2.2) \x\ = sup | Aj|. 

i 

In the case of Sym(r, R), (|2.ip corresponds to the polar decomposition 
(diagonalization) of symmetric matrices. The functions det and tr are 
respectively the usual determinant Det and trace Tr of matrices. 



3. Hilbert's metric on symmetric cones 

If we consider the cone f^sym of real symmetric positive definite r xr 
matrices, then one can easily express the Hilbert projective metric (II. lj) 
in terms of eigenvalues of elements of fisym- Indeed, if A and B are in 
fi Sym , then 

M(A, B) := inf{A \ XB - A 4 0} = max [ Ax \ x \ , 

H=i \Bx \x ) 
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and 



m{A, B) := sup{A | XB - A ^ 0} 



(Ax\x) 

= mm 

IMI=i [Bx \x) 



which are respectively the greatest and the least eigenvalue of B~ l A. 
Observe that eigenvalues of the matrix B~ l A are the same of matrix 
B~zAB~2 = P(B~z)A (see Example 12. II for this notation). 

More generally, for symmetric cones, Hilbert's projective metric can 
be also formulated in terms of extremal eigenvalues : let x and y be in 
Q and let Aa/(x,?/) > and X m (x,y) > denote the greatest and the 
least eigenvalue of the element P(y~^)x G Q. Then one can prove the 
following (see [7| Thoerem 4.2]) 

Proposition 3.1. We have 



From this characterization one can easily prove the following 

Lemma 3.2. If x, y 6 and a, (3 > 0, then 

(1) X M (ax + fly, y) = a\ M (x, y) + (3, 

(2) \ m (ax + f3y, y) = a\ m (x, y) + (3, 

(3) \ M (x,y)\ m (y,x) = 1. 

Proposition 3.3. If x, y e Q, then the Hilbert metric of x and y is 



(3.1) 



^M{x,y) 




and 



(3.2) 



Kn{x,y) 





Proof. According to Proposition 13. 11 we have 
M(x,y) = inf{A | x ^ Xy} 



= inf{A \ Xy — x E Q} 
= m^{^\ceJ(V)} 
= X M (x,y), 
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Applying LemmaHOl we obtain easily the second equality of the propo- 
sition. □ 

Remark 3.4. It follows from Lemma f3. 21 and Proposition 13.31 that d is 
constant on rays, 

(3.4) d(Xx, /xy) = d(x,y) for A, \i > 0. 

In particular, if \x > 0, then the map ^ : x i — > fix is an isomery of d. 

4. Contractions of Hilbert's projective metric 

Let x = X/*=i ^j c j be an element of Q. If p G R, we write x p for 
Sj=o ^j c j anc ^ we state the following contraction result : 

Proposition 4.1. Let p G R such that \p\ < 1. Then the map 

(4.1) cl> p : Q — > fi, xhj p 
is a contraction of d, that is, for any x, y G f2, 

(4.2) g?(x^) < \p\d(x,y). 

Proof. Let x, y G Q. Then Am(^ _1 , Z/ -1 ) = 1/A m (x, y) and A m (x _1 , y _1 ) 
1/Xm(x, y). Hence the map 

(4.3) l: ll^fl:iHi 4 

is an isometry of d. It suffices then to prove ()4.2j) for < p < 
1. We know that x ^ \M{x,y)y, then by the Loewner-Heinz in- 
equality, see p, Corollary 9], we have x p ^ XM(x,y) p y p . Therefore 
Am(^ p ,2/ p ) < Am(^c, y) p - Using the same arguments we prove that 
Xm(x p ,y p ) > Xm(x,y) p . Thus d{x p , y p ) < pd{x, y). □ 

5. Completeness 

Proposition 5.1. (Q, d) is a pseudo-metric space. In other words, for 
any x, y, z G Q, the following holds : 

(a) d(x,y)>0 

(b) d(x,y) = d{y,x) 

(c) d(x,z) < d(x,y) + d(y,z) 

(d) d(x, y) = <^> 3A > : x = \y 

Proof. Let x, y and z be in f2. The first propriety of the proposition is 
obvious, since \M{x,y) > X m (x,y). 
(b) From Lemma f3. 21 we obtain 

= logX M (x,y)X M (y,x) 
= d(y,x). 
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(c) Let c be a primitive idempotent of V, then according to Proposition 
O < gj| < X M (x,y) and < g| < \ M (y,z). Hence gjf = 

§§{5| < A m (^,2/)Am(2/,^) and 

(5.1) Xm(x, z) < Xm{x, y)X M (y, z). 
Similarly we prove that 

(5.2) X m (x,z) > X m (x,y)X m (y,z). 

It follows then from (|5.1|l and ()5.2|) that d(x, z) < d(x, y) + d(y, z). 

(d) If d(x,y) = 0, then Am(,^,2/) = X m (x,y) := X and all the eigen- 
values of P(y~2^x are equal to A. Therefore, P(y~z)x = Xe and 
x = XP(yz)e = Xy. Conversely, if x = Xy where A > 0, then 
P(y~2)x = Xe and X M (x,y) = X m (x,y) = X. □ 



Let S(V) be the unite sphere in V, S(V) = {x £ V \ \x\ = 1} with 
respect to the spectral norm introduced in (|2.2j) . 

Lemma 5.2. Let x and y be in Q(l S(V). Then 

(a) \x-y\< e d{x ^ - 1. 

(b) If \x — y\ < X m (y) then \x — y\ > X m (y) tanh{|d(x, where 
X m (y) is the least eigenvalue ofy. 

Proof. Let x,y £ On S(V) and let c be a primitive idempotent of V, 
then by Proposition (jH.lj) we have 

(5.3) A m (x,y)(y|c) < (x\c) < X M (x,y)(y\c). 

But |x| = 1 and \y\ = 1, then there exists c\,c% £ J{V) such that 
(x|ci) = 1 and (y\c 2 ) = 1, thus 

(5.4) X m (x,y) < 1 < X M (x,y). 

Let c be a primitive idempotent of V such that |rr — y\ = (x — y\c ), 
then using ()5.4|) we have 

\x-y\ = (x\c ) - (y\c ) 

= (S-l}(^o) 

< X M (x,y) - 1 

< X M (x,y) - X m (x,y) 

< {X M (x,y) -X m {x,y)}j^Tyj 
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Moreover 

X M {x,y) = max — — 
ceJ(V) (y\c) 

(x-y\c) 
= max { — — — hi} 

c€j(v) 1 (y\c) 

(5.5) < \x-y\/X m (y) + l. 

Now, if \x — y\ < X m (y), then 1 — > ®- Therefor, from (|5.5j) 

and Lemma l'A.'2\ we have (l — jp~7~l ) V x. Hence 

(5.6) \ m {x,y)>l-^L. 

Finally, it follows from (|5.5|) and (|5.6|) that \x—y\ > X m (y) tanh{lc?(x, y)}. 

□ 

Proposition 5.3. (Q n S(V),d) is complete metric space. 

Proof. It is clear that if x, y G (f2nS'(V A ), d) such that y) = 0, then 
x = y. Thus (Q D S'(V), d) is a metric space. 

Let (a;jfc)fe be a Cauchy sequence in (QnS(V), d). Then from Lemma l5~2"l 
(xfc)fc is a Cauchy sequence in (V, | • |) and then converges to an element 
iGfln 5'(V A ). In order to prove that x is an element of SI Pi S(V), we 
observe, using (|5.2|) that A m (xj) = A m (rr i3 e) > X m (xi,Xj)X m (x j: e) = 
X m (xj)XM(xi, Xj)e~ d( - Xl ' Xj \ Therefor, from ()5.4j) . we obtain A m (xj) > 
X m (xj)e~ d ( Xi > Xj \ Since (xfc)& converges in the spectral norm to x, (X m (x k )) k 
converges to X m (x). It follows that for a fixed large j, X m (x) = 
lim^oo X m (xi) > aX m (xj) > 0, where a > 0. Hence x G Q n S'(V). 
Finally, it is easy to prove using the second proposition of Lemma 15.21 
that (xk) converges to x in (O fl -S'(V), <i). □ 



6. Application 

If V is the Jordan algebra Sym(r, R), then Q is the symmetric cone 
of positive definite matrices, and G(Q) coincides with the linear group 
GL(r, R). Recall that in this case G(Q) acts on Q by : t ■ a = t 'at. 
In |3] , Bushell prove that for any t G GL (r, R) and any integer > 1 , 
there exists a unique real positive definite matrix a such that if at = a 2 . 
For a general symmetric cone we have the following 

Theorem 6.1. Let g G G(Q) and p G R st/c/i t/iat \p\ > 1. T/ien t/iere 
exists a unique element a in Q such that g(a) = a p . 
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Proof. Since i : x i— > x 1 is an isometry of d, see f!4.3|l . we can assume 

that p > 1. For s G fl S'(V) we put = — ^^(7(2;)?. Then F 

Is0) p l 

maps fl S(V) into fl S(V). Recall that Hilbert's projective metric 
is constant on rays, see (|3.4|) . then we have for x, y G Q fl S(V), 

d(F(x), F{y)) = d((wi/ p o g)(x), o g)(y), 

where (jJ\i v is the power map x t— > x?. We have already noticed in (j4.1|) 
that ui/p is a contraction of the metric d, hence 

d(F(x),F(y))<-d(g(x),g(y)). 
V 

Moreover, the eigenvalues of P(y~ l l 2 )x are the unique solutions of the 
characteristic equation det(Ay — x) =0. But, using [5, Proposition 
III.4.3], 

&et{\g{y) - g{x)) = det{g{\y - x)) = Vet{g) n ' r det{\y - x). 

Thus the eigenvalues P(y)~ 1 ^ 2 x and P(g(y))~ 1 ^ 2 g(x) are the same and 
g is also an isometry of d. Finally all this implies that F is a contraction 
of the Hilbert metric d, 

d(F(x),F(y))<-d(x,y). 

V 

Therefore by the Banach contraction mapping theorem and Proposi- 
tion 15.31 the map F has a unique fixed point, let say u. Then u p = 
\ g (u)i/p\p 9( u ) an d the element a := \g(u) l l p \ p u G Vt is the unique solu- 
tion of the equation g(x) = x p . □ 

Corollary 6.2. Let h G G(Q) and pel such that \p\ > I. Then 
there exists a unique element a of Q such that h(a p ) = a. 

Final remark. The Hilbert original definition of the projective metric 
involved the logarithm of the cross-ratio of for points in the interior of 
a convex cone in IR n . It would be interesting to define the Hilbert 
projective metric on symmetric cones using the generalized cross-ratio 
introduced in [2], see also jTj. 
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